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The estimated brightness of the Cygnus loop supernova remnant (SNR) at 2720, 1420, 820, 408 
and 34.5 MHz are presented. The observations of the continuum radio emission are used to calculate 
the mean brightness temperatures and surface brightnesses of this loop at the five frequencies in wide 
spectral range, using the method we have previously developed for large radio loops. The spectrum 
for mean temperatures versus frequency between the five frequencies is estimated and the spectral 
index of Cygnus loop is also obtained. Also, from our results can be concluded that Cygnus loop 
evolves in the low density environment and the initial energy of supernova explosion was relatively 
low. The obtained results confirm non-thermal origin of the Cygnus radio loop and show that our 
method is applicable to almost all remnants. 
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I. INTRODUCTION 

In our investigation, we are supposing the radio loops 
to be evolved supernova remnants (SNRs). Super- 
nova remnants represent shells of approximately spher- 
ical shape which are spreading, and depending of the 
interstellar matter density, they change their shape. We 
can see only parts of these spherical contours which rep- 
resent radio spurs, e.g. areas of more intense emission in 
the sky, spur-like and of huge dimensions. Radio loops 
consist of more spurs lying approximately in the same 
small circle of the celestial sphere. Their material prob- 
ably expands inside of bubbles of low density. Bubbles 
are made by former SNR explosions or by strong stellar 
winds (McKee & Ostriker [l(, Salter Q and references 
therein). The radio emission from SNRs is generally un- 
derstood to be synchrotron emission from the relativistic 
electrons moving in the magnetic field. 

A star in the constellation of Cygnus exploded and its 
remnant is Cygnus loop. It is classified as a middle-aged 
SNR of type S, located below (but near the plane of) 
the Galactic equator, and less than 1 kpc away from us. 
It is listed in Green's catalogue of SNRs as G74.0-8.5 
(Green 0,1, H]). As shown in Green Q, this remnant 
has been decelerated considerably by its interaction with 
the surrounding interstellar medium. 

Leahy, Roger, & Ballantyne Q presented high resolu- 
tion 1420 MHz total intensity and polarization maps of 
the Cygnus loop and derived rotation measures. Aschen- 
bach & Leahy Q made a comparison between the radio 
and X-ray emission from supernova remnants. The X-ray 
emission traces the hot gas behind the shock front, while 
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the radio emission comes from the relativistic electron 
population emitting synchrotron radiation in the ambi- 
ent magnetic field. They found that there are significant 
differences in the distribution of X-ray and radio bright- 
ness in the Cygnus loop. Also, significant temperature 
variations, seen along the rim and in the bright filaments 
interior to the rim, pointed out that the Cygnus loop was 
caused by a supernova exploding in a cavity. 

Uyaniker et al. Q suggested that Cygnus loop may 
consist of two overlapping remnants. Substantial differ- 
ences between the northern and southern part as well as 
in their emission characteristics have been observed (Pat- 
naude et al. [13], Uyaniker et al. [11])- Some authors 
proposed that the Cygnus loop consists of two likely in- 
teracting SNRs: G74.3-8.4 and G72.9-9.0 (Uyaniker et 
al. @, Leahy [12]). In Sun et al. 0, from 4800 MHz 
observations, it is explained that the polarization maps 
(the difference in the polarization characteristic between 
the northern and southern part) support previous ideas 
that the Cygnus loop may consist of two SNRs. Be- 
sides, several compact radio sources are located within 
the boundary of the remnant. The main characteristics 
of the Cygnus loop in different spectral bands are: (a) 
in radio band: it is shell, brightest to the north-east, 
with fainter breakout region to south, with spectral vari- 
ations, (b) optical: large filamentary loop, brightest to 
the north-east, not well defined to the south and west, 
(c) X-ray: shell in soft X-rays (see e.g. in Green Q). 

Some radio maps which include this remnant are the 
following: the map by Kundu & Becker [ijj at 4940 
MHz, Keen et al. [ll at 2695 MHz, Uyaniker et al. 
Tl| at 2675 MHz, Leahyjloger, & Ballantyne at 1420 
MHz, Dickel & Willis O at 610 MHz, Green [| at 408 
MHz. Observations of the continuum radio emission at 
2720 (Reif et al. [13), 1420 (Reich & Reich [13), 820 
(Berkhuijsen 408 (Haslam et al. (U) and 34.5 

MHz (Dwarakanath & Udaya Shankar [li[), we found in 
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electronic form and used them in this paper. 

In our previous paper (Borka Jovanovic & Urosevic 
[22j) we only calculated the temperatures and bright- 
nesses at the three frequencies. In this paper we expand 
the scope of our investigation to: brightness tempera- 
tures and surface brightnesses at the five frequencies, as 
well as spectrum, T— T graphs, the radio spectral indices, 
estimation of environment density and the initial energy 
of supernova explosion, and the flux density spectrum. 

In our calculations, we used brightness temperatures 
over the whole area of the loop, so the mean temperature 
that we estimated refers to northern and southern parts 
together. In this research the average brightness tem- 
peratures and surface brightnesses of the Cygnus radio 
loop are calculated at the five frequencies: 2720, 1420, 
820, 408 and 34.5 MHz. Then we study how these re- 
sults are getting along with previous results (Uyaniker et 
al. [U, Reich, Zhang & Fiirst [23], Roger et al. [H|) 
and with current theories of SNR evolution. These the- 
ories predict that loops are non-thermal sources which 
are spreading inside of the hot and low density bubbles 
made by former supernova explosions or by strong stellar 
winds (see McKee & Ostriker [l|, Salter [2| and refer- 
ences therein). 

Our aim is also to apply method for determination of 
the brightness temperature given in article Borka [25[ 
which is developed for large radio loops, and to show 
that it is rather efficient in the case of much smaller ra- 
dio loops, e.g. Cygnus loop, as well as to check how the 
results obtained using this method are getting along with 
the method of T — T graphics and results obtained with 
other methods. Our method is quite simple because we 
are using brightness temperature isolines to define bor- 
ders of the Cygnus loop at the wide range of frequencies 
{Vmax/Vmin = 2720 MHz / 34.5 MHz w 79). Other au- 
thors are using different squared or rectangular areas to 
determine area of the loop and calculate spectral indices, 
brightness temperature and the flux of the loop (Uyaniker 
et al. 11 1, page 917 and Leahy & Roger [26], page 786). 
Also, we calculated flux from Cygnus loop and compare 
our results with results of other authors in different fre- 
quency ranges. 



II. ANALYSIS 
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Data 

We used the following radio-continuum surveys as the 
basic source of data in this paper: at 2720 MHz (Reif 
et al. [13), 1420 MHz (Reich & Reich 0), 820 MHz 
(Berkhuijsen Qjj), 408 MHz (Haslam et al. [H) and 
34.5 MHz (Dwarakanath & Udaya Shankar [13). These 
surveys are available in electronic form, as the "Flexible 
Image Transport System" (FITS) data format, at the 
site: |http : / / www . mpif r-bonn . mpg . de /survey . html| 
This online Survey Sampler of the " Max-Planck-Institut 
fur Radioastronomie" (MPIfR) near Bonn, Germany, al- 
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FIG. 1: Top: the 2720 MHz map of a region in Cygnus, in new 
Galactic coordinates (/, 6), showing contours of brightness 
temperature. This radio loop has position: I — [76.5°, 71.5°]; 
b = [-10.5°, -7°]. The HPBW (Half Power Beam Width) for 
this frequency is 0°.35. Eleven contours plotted represent the 
temperatures T m i n and T max from Tableland nine contours 
in between. The contours are plotted every 0.265 K, starting 
from the lowest temperature of 0.395 K up to 0.66 K. The 
corresponding temperature scale is given (in K). Bottom: the 
2720 MHz area map of Cygnus. 



lows users to pick a region of the sky and obtain images 
and data at different frequencies. The 2720-MHz Stock- 
ert survey has the angular resolution of 0°.35, 1420-MHz 
Stockert survey (Reich & Reich [lil)0 o .59, the 820- 
MHz Dwingeloo survey (Berkhuijsenjll) 1°.2, the 408- 
MHz all-sky survey (Haslam et al. [20]) 0°.85 and the 
34.5-MHz Gauribidanur survey (Dwarakanath & Udaya 
Shankar [2l|) 0°.7. The corresponding observations are 
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FIG. 2: The same as Fig. 1, but for 1420 MHz. The HPBW 
for this frequency is 0°.59. The contours are plotted every 
0.07 K, starting from the lowest temperature of 4.2 K up to 
4.9 K. 




FIG. 3: The same as Fig. 1, but for 820 MHz. The HPBW 
is 1°.2. The contours are plotted every 0.39 K, starting from 
the lowest temperature of 10.1 K up to 14 K. 



given at the following rates (measured data) for both I 
and b: ^ at 2720 MHz, ^ at 1420 MHz, ^ at 820 

MHz, 2£ at 408 MHz and ^ at 34.5 MHz. The effec- 
tive sensitivities are about 5 mK Tt (T& is for an average 
brightness temperature), 50 mK, 0.20 K, 1.0 K and about 
700 K, respectively. 

FITS data format stores the multidimensional arrays 
and 2-dimensional tables containing rows and columns 
of scientific data sets. We extracted observed brightness 
temperatures from this data format into ASCII data files, 
and afterwards using our programs in C and FORTRAN, 
we obtained results presented in this paper. 



Method 

As there is great influence of background radiation over 
the Cygnus loop area, it is very difficult to determine 
its area precisely. The area of this loop is enclosed with 
brightness temperature contours. The maps of a region in 
Cygnus, in new Galactic coordinates (I, b), with contours 
of the brightness temperatures It, are plotted in Figs. [1|- 
[5] In these figures, among all the brightness temperature 
contours, the most important are the outer contour and 
inner one which represent the loop borders. The outer 
one (which corresponds to the minimum temperature of 
the loop) separates loop and background, while the inner 




T(K) 




FIG. 4: The same as Fig. 1, but for 408 MHz. The HPBW 
is 0°.85. The contours are plotted every 3.1 K, starting from 
the lowest temperature of 41 K up to 72 K. 



one (corresponding to the maximum temperature) sepa- 
rates loop and some superposed source. The interval of 
Galactic longitude of this loop is I = [76.5°, 71.5°], and 
of latitude is b — [-10.5°, -7°]. We used the same method 
of calculation as given in Borka [25j for Galactic radio 
loops I- VI, Borka, Milogradov- Turin, & Urosevic 27 1 
for Loops V and VI and in Borka Jovanovic & Urosevic 



[28j | for Monoceros loop. Our aim is to apply method, 
which we developed for main Galactic loops I- VI and de- 
scribed in Borka [25| and Borka, Milogradov- Turin, & 
Urosevic [13] , to smaller remnants and to show that our 
method of calculation is applicable to almost all rem- 
nants. To confirm that our method is good, we also gave 
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FIG. 5: The same as Fig. 1, but for 34.5 MHz. The HPBW 
is 0°.7. The contours are plotted every 3550 K, starting from 
the lowest temperature of 34500 K up to 70000 K. 



the area maps of Cygnus at the five frequencies (Figs. 
HHS]) and some examples of the temperature profiles at 
the 2720 MHz (Fig. E]). In Fig. [6] it should be noticed 
that the bottom panel shows only background radiation 
(not including loop), where for b = -10° the brightness 
temperature is not higher than 0.38 K. It is in agreement 
with temperature intervals given in the first row of the 
Table U where for the minimum temperature of the loop 
with background, at 2720 MHz, it is put 0.395 K. 
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FIG. 6: Temperature profiles at 2720 MHz for Galactic longitude from 76°. 5 to 71°. 5 and for the following values of Galactic 
latitude: b = -7°. 5 (top left), b = -8° (top right), b = -8°. 5 (middle left), b = -9° (middle right) and b = -10° (bottom). 
Temperatures are given in K, and Galactic longitudes in degrees. Notice that the bottom figure shows the temperature profile 
at latitude which lies outside (but close enough to) the loop and from that profile we can see what is the highest temperature 
of the background at b — -10°, at 2720 MHz (not of the loop with background). 
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The mean temperatures and surface brightnesses of 
this radio loop are computed using data taken from radio- 
continuum surveys at 2720, 1420, 820, 408 and 34.5 MHz. 
We have subtracted the background radiation, in order 
to derive the mean brightness temperature of the SNR 
alone. First, the temperature of the loop plus background 
was determined. For every bin in survey we have some 
value for temperature T sum . We take all bins within the 
loop border and obtained average temperature < T sum > 
of the loop plus background. Then, the background alone 
near the loop was estimated. We determine temperature 
of the background near the outer border of the loop (aver- 
age value from all data near the outer border of the loop) 
and that is the background temperature Tback- And fi- 
nally, the difference of these values is calculated to obtain 
the average temperature of the loop at each frequency: 
average temperature of the loop is < T surn > —Tback- 
The areas over which an average brightness temperature 
is determined at each of the five frequencies are taken 
to be as similar as possible within the limits of mea- 
surement accuracy. T m ; n from Table U means the lower 
temperature limit between the background and the loop, 
and T max means the upper temperature limit of the loop. 
So we used all measured values between T min and T max , 
inside the corresponding regions of I and 6, to calcu- 
late the brightness temperature of a loop including the 
background. The mean brightness temperature for the 
loop is found by subtracting the mean value of back- 
ground brightness temperature from the mean value of 
the brightness temperature over the area of the loop. 

After deriving the mean brightness temperatures T\, tl/ , 
we have converted these values into surface brightness £„ 
by: 



2kv z /c') T h>v , 



(1) 



where k is Boltzmann constant and c the speed of light. 
Results are given in Table fl] 



are defined to separate the spur and its background. For 
the method of calculation see also Borka [25]], Borka, 
Milogradov- Turin, & Urosevic [13] and Borka Jovanovic 
& Urosevic |28[. As mentioned in these papers, if the 



value of T m j n is changed by a small amount, the bright- 
ness contours become significantly different. If T m ; n is 
too small, the area of the spur becomes confused with 
the background and it becomes obvious that the border 
has been incorrectly chosen. 

The results are given in Tableslll iHland lHll T m j n , given 
in the second column of Table U is the lower temperature 
limit, while T max is the upper temperature of the loop 
and it is also upper temperature limit (because there are 
no other superposed sources). These temperature limits 
enable us to distinguish the loop from background. Then 
we derived the surface brightnesses using equation ([I} for 
each frequency. 

The values for brightnesses in 10~ 22 W/(m 2 Hz Sr) can 
be compared with results for flux densities in Jy. The 
flux densities can be transformed into brightnesses or vice 
versa, taking into account that the frequencies have to be 
the same. Knowing the loop size f2, the flux densities S v 
given in Jy can be transformed to brightnesses given in 
10~ 22 W/(m 2 HzSr) by: 



(2) 



By use of the spectral indices, the brightnesses can be 
reduced to 1000 MHz according to relation: 



Siooo = £„(1000/i/) 



(2-/3) 



(3) 



where the temperature spectral index (3 — a + 2, and a 
is the spectral index defined by S v oc v~ a . 



Spectrum 



III. RESULTS 

The radio continuum maps are used for determining 
the area of the Cygnus loop and for deriving brightness 
temperatures over it. At each of the five frequencies, the 
areas are determined to be as similar as possible within 
the limits of measurement accuracy. There are still some 
differences between these areas and we think that the ma- 
jor causes of differing borders between the five frequencies 
are small random and systematic errors in the data. The 
surface brightnesses of SNRs must be above the sensitiv- 
ity limit of the observations and must be clearly distin- 
guishable from the Galactic background emission (Green 
[29j|). As it is very difficult to resolve the fainter parts 
of the loop from the background, they are not taken into 
account. For evaluation brightness temperatures over the 
area of the loop we had to take into account background 
radiation (see Webster [3(| )■ Borders enclosing the spurs 



The spectrum was generated using mean temperatures 
at five different frequencies. Best-fit straight line spec- 
trum enables calculation of spectral index as negative 
value of the line's direction coefficient. In Figure [7] we 
give spectrum for the three middle frequencies: 1420, 
820 and 408 MHz because their spectrum is very well 
fitted with the straight line (see Figure [7]) . Frequencies 
2720 and 34.5 MHz, lies on very high and on very low 
ends of the spectrum, as presented in Fig. [5J If for cal- 
culation of spectral index we take only three frequencies 
1420, 820 and 408 MHz, the result would be /3 3 = 2.76 
± 0.03. All five frequencies 2720, 1420, 820, 408 and 
34.5 MHz, also from linear fit, give (3 5 = 2.66 ± 0.09. It 
is steeper spectral index in comparison to average value 
for Galactic SNRs (3 = 2.5 (in Green's catalogue (Green 
Q) the adopted value for spectral index is not given be- 
cause of its variability). Cygnus loop is relatively old 
SNR. The shock wave has to be weak for evolved SNRs. 
The steeper spectral indices of SNRs should be expected 
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TABLE I: Temperatures and brightnesses of Cygnus radio loop at 2720, 1420, 820, 408 and 34.5 MHz 



Frequency 


Temperature limits 


Temperature 


Brightness 


(MHz) 


Tmin, ^max (-^) 


(K) 


(10~ 22 W/(m 2 Hz Sr)) 


2720 


0.395, 0.66 


0.160 ± 0.005 


3.64 ± 0.12 


1420 


4.2, 4.9 


0.49 ± 0.05 


3.04 ± 0.30 


820 


10.1, 14.0 


2.30 ± 0.20 


4.75 ± 0.40 


408 


41, 72 


15.3 ± 1.0 


7.83 ± 0.50 


34.5 


34500, 70000 


13960 ± 700 


50.97 ± 2.56 



for the smaller shock wave velocities. It is result ob- 
tained from the diffuse particle acceleration theory (Bell 
[3lT Hl|)- Additionally, Cygnus loop probably expands 
in low density environment. It looks like large Galac- 
tic radio loops - evolved SNRs with the steep spectral 
indices, immersed in the low density environment (see 
Borka 



251 ] . Borka, Milogradov- Turin, & Urosevic 



and references therein). 




logv 

FIG. 7: Cygnus loop spectrum: temperature versus frequency, 
for three measurements - at 408, 820 and 1420 MHz. 

Obtained values ^3 = 2.76 ± 0.03 (from three fre- 
quencies), /?5 = 2.66 ± 0.09 (from all five frequencies), 
are greater than 2.2 and confirm non-thermal origin of 
Cygnus loop emission. From Fig. |S]it can be noticed that 
linear fit is quite satisfactory. The value for the bright- 
ness temperature spectral index of the Cygnus loop is 
rather steep. This is at the high end of the spectral index 
distribution for SNRs as suggested in Clark & Caswell 

M. 



T -T plot 

The measured data have different resolutions for 
different frequencies (see §2.1), and therefore in order 
to obtain T — T plots the data were retabulated so the 
higher resolution maps are convolved to the resolution 
of the lowest resolution map. In that way we convolved 
data at 2720, 1420, 408 and 34.5 MHz to 0°.5 x 0°.5 




log V 



FIG. 8: Cygnus loop spectrum: temperature versus frequency, 
for five measurements - at 34.5, 408, 820, 1420 and 2720 MHz. 



resolution, which is the sampling rate of the 820 MHz 
survey. These retabulated data are presented in Figure 
OH for the following frequencies: 2720 MHz (top left), 
1420 MHz (top right), 820 MHz (middle left), 408 MHz 
(middle right) and 34.5 MHz (bottom). Then, for each 
frequency pair we used only the common points (with 
the same (I, b)) which belong to the loop area at both 
frequencies. In that way we reduced loop area to the 
same area for different frequencies. The obtained T — T 
plots for five pairs of frequencies enabled calculating the 
spectral indices. We calculated two (3 values for each of 
these pairs: between 2720-1420, 2720-820, 2720-408, 
2720-34.5, 1420-820, 1420-408, 1420-34.5, 820-408, 
820-34.5 and 408-34.5 MHz and presented it in TableHU 
For each of the ten frequency pairs, by interchanging the 
dependent and independent variables we have obtained 
two /3 values for each pair and the mean value of these 
fit results is adopted as the radio spectral index, as 
suggested in Uyaniker et al. [ll|. Regarding only three 
frequencies (because their spectrum lies on straight line, 
see Figure [7} 1420, 820 and 408 MHz, the average value 
of spectral index from T - T is < fi TT > 3 = 2.78 ± 
0.41. Taking into account all five frequencies, we get 
< Ptt >5 = 2.63 ± 0.30. It can be noticed that this 
value agrees well with the corresponding value obtained 
from spectrum, as expected (see Uyaniker et al. llj). 




FIG. 9: The data retabulated to 0°.5 x 0°.5 resolution, for the following frequencies: 2720 MHz (top left), 1420 MHz (top right), 
820 MHz (middle left), 408 MHz (middle right) and 34.5 MHz (bottom). The HPBWs for these frequencies are 0°.35, 0°.59, 
1°.2, 0°.85 and 0°.7, respectively. The gray scales of temperatures are given also. 
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TABLE II: Spectral indices for Cygnus loop from T — T plots, between 2720, 1420, 820, 408 and 34.5 MHz 



Frequency (MHz) 


2720 


1420 


820 


408 


34.5 


2720 


/ 


1.52 ± 0.35 


1.96 ± 0.53 


2.35 ± 0.26 


2.57 ± 0.23 


1420 


2.22 ± 0.76 


/ 


2.27 ± 0.64 


2.60 ± 0.28 


2.67 ± 0.25 


820 


2.64 ± 0.61 


3.10 ± 0.64 


/ 


2.82 ± 0.31 


2.82 ± 0.10 


408 


2.67 ± 0.26 


2.79 ± 0.28 


3.10 ± 0.31 


/ 


2.78 ± 0.15 


34.5 


2.86 ± 0.23 


2.95 ± 0.25 


2.92 ± 0.10 


2.94 ± 0.15 


/ 



Then we calculated mean value of spectral index: re- 
garding spectrum and T — T graphs. Between 1420, 820 
and 408 MHz we obtained < /3 > 3 = 2.77 ± 0.22, and 
between 2720, 1420, 820, 408 and 34.5 MHz we obtained 
< p > 5 = 2.64 ± 0.20. 



IV. DISCUSSION 

Spectral index variations with position of more spatial 
features within the Cygnus loop can be found in paper 
Leahy (34[. There were studied radio spectral indices 
by T - T plot method between 2695, 1420 and 408 MHz 
and found that the bright radio filaments all show nega- 
tive curvature (steeper at higher frequency), and regions 
dominated by diffuse emission show positive curvature 
(flatter at higher frequency). 

It can be noticed that mean value of spectral index for 
Cygnus loop is little higher then the corresponding value 
obtained from articles (see Uyaniker et al. [ll[ or Leahy 
& Roger (26|). Reason for this is different areas used for 
data of Cygnus loop in these papers. In both papers they 



11] 



used square areas (see page 917 from Uyaniker et al 
and page 786 from Leahy & Roger [26J ) . In these square 
areas they take into account parts of Cygnus loop and 
parts of background radiation near Cygnus loop (these 
parts we do not take into account). We take in account 
only higher intensity regions from these areas (see Figures 

HHS1). 



TABLE III: Brightness of Cygnus radio loop reduced to 1000 
MHz, using spectral index derived in this paper: /?5 = 2.66 
± 0.09 (from the spectrum of all five frequencies 2720, 1420, 
820, 408 and 34.5 MHz) 



Frequency 


Brig 


htness at 1000 MHz 


(MHz) 


(10 


" 22 W/(m 2 Hz Sr)) 


2720 




7.05 ± 0.86 


1420 




3.84 ± 0.52 


820 




4.17 ± 0.29 


408 




4.33 ± 0.07 


34.5 




5.52 ± 1.40 



In order to compare our values for brightnesses with re- 
sults for fluxes given in other papers, we transform their 
calculated fluxes into brightnesses at 1000 MHz. Know- 
ing the loop size fl we reduce the flux densities given in Jy 
to brightnesses given in 10 -22 W/ (m 2 Hz Sr) by equation 
@ and then by use of the spectral indices, the bright- 



nesses are extrapolated to 1000 MHz according to rela- 
tion ©. 

In our previous paper (Borka Jovanovic & Urosevic 
[22|| ) we transformed the flux densities for Cygnus Loop 
given in Roger et al. El at 22 MHz (S v = 1378 Jy) and 
Reich, Zhang & Fiirst [23| at 863 MHz (S v = 184 Jy) into 
E1000 using loop size ft — 240' x 170' from Reich, Zhang 
&; Fiirst [23[ and spectral index /3 = 2.49 from Trushkin 
[35l | . We obtained that values from (Borka Jovanovic & 
Urosevic [22]) agree with previous data. But now in 
this paper, we use our derived value for spectral index 
(3$ = 2.66±0.09, and obtain even better agreement. From 
flux densities given in mentioned papers, we calculated 
the following values for radiation intensities: S1000 = 
3.22 x 10~ 22 W/(m 2 HzSr) from flux given in Roger et 
al. [H and Ei 00 o = 4.84 x 10~ 22 W/(m 2 Hz Sr) from 
flux given in Reich, Zhang & Fiirst [23] ]. 

Using our calculated (3$ and applying relation ([3]) 
to our calculated values of E„ at 2720, 1420, 820, 
408 and 34.5 MHz, for S1000 wc obtain values given 
in Table IIIII Absolute errors for brightnesses at 
1000 MHz are calculated in this way: ASiqoo ~~ 
(Eiooo/S„i) (AE„! + E„i • In (i^/IOOO) • A/3), where v Y 
takes values 2720, 1420, 820, 408 and 34.5 MHz. 



TABLE IV: The Cygnus loop areas Q and its angular radii 
at the five frequencies 



Frequency (MHz) 


Q ((T) 


en 


2720 


7.52 


1.547 


1420 


8.64 


1.658 


820 


9.46 


1.735 


408 


9.88 


1.773 


34.5 


8.44 


1.639 



Information about variations of brightness and spec- 
trum over radio image is presented by Leahy & Roger 
[26l | and by Uyaniker et al. [9]. The average values re- 
duce the information imprinted in variations of bright- 
ness and spectrum over the radio image but give us 
other interesting information, like explosion energy, sur- 
face brightness and distance to the Cygnus radio loop. 

With regard to the loop's borders in Figs. HHU we 
derived loop area at each frequency, as well as its cor- 
responding angular radii (see Table IIV[) . An angular 
radius is obtained in this way: our derived loop area, 
when approximated with the circle of the same area, gives 
possibility of determining angular radius as 9 = y^T/V. 
These areas can be compared to the areas calculated by 
other authors, e.g. 9(°) 2 .80 at 1420 MHz (Leahy [Hj]) 
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FIG. 10: The surface brightness to diameter diagram from 
Berezhko & Volk (2004), with value for Cygnus loop added. 
Three different densities for the interstellar matter (ISM) (A*h 
= 3, 0.3 and 0.003 cm -3 ) are presented, plus three values for 
the (SN) explosion energy (E S n = 0.25, 1 and 3 xlO 51 erg). 



and 11(°) 2 .33 at 863 MHz (Reich, Zhang & Furst H). 
When their areas are recalculated into angular radii in 
the way we described, the result are: 1° .77 at 1420 MHz 
(Leahy El) and 1°.90 at 863 MHz (Reich, Zhang & 
Furst |23||). It can be noticed that we obtained some- 
what smaller values, but it has to be taken into account 
that other authors estimated only the rectangular map 
size while we estimated the loop size (inside its contour 
borders) exactly. 

The relation between surface brightness (E) and diam- 
eter (D) for supernova remnants (SNRs) - so-called the 
E— D relation - is appropriate for description of the radio 
brightness evolution of these sources. The empirical rela- 
tions should be used with caution because of the limited 
usefulness due to selection effects (see Green [j| , Urosevic 
et al. [IB] and references therein). The updated theo- 
retical relations were derived by Duric & Seaquist [37j . 
based on the Bell's ([3l|, HH) diffuse shock particle ac- 
celeration theory, and by Berezhko and Volk (Berezhko 
& Volk (38|) based on the non-linear kinetic theory of 
the diffuse shock acceleration mechanism. The E — D 



diagram at 1 GHz with the theoretically derived evolu- 
tionary tracks taken from Berezhko & Volk [38| with our 
value for the Cygnus loop superposed, is shown in Fig. 
ITU1 The E — D diagram at 1 GHz taken from Berezhko 
& Volk [38[ with the derived value for the Cygnus loop 
superposed, is shown in Fig. [IJJ With aim to superpose 
Cygnus position, we used the following results of our cal- 
culations: the mean value of brightnesses at 1 GHz and 
the diameter, calculated as D = 2rsin(9, where 9 = 1°.67 
being mean angular radius for all five frequencies. The 
distance r — 0.44 kpc is taken from Green's catalogue of 
SNRs (Green @). So we added this point to the dia- 
gram: (D, E) = (25.7 pc, 4.98 x 10~ 22 W/(m 2 Hz Sr)). 
From its position in this diagram it can be concluded 
that Cygnus loop evolves in the low density environment 
and the initial energy of supernova (SN) explosion was 
relatively low (see Fig. [T0|) . 



Flux density spectrum 

On the basis of values for brightnesses given in the 
fourth column of Table fl] and our calculated values for 
loop size f2 (Table IIVI) . we derived flux values in Jy. 
The calculated flux densities we give in Table IVl As the 
Cygnus loop is well studied SNR, it is possible to per- 
form a multi- frequency spectral study. In the Fig. [11] we 
present a summary of some results obtained from other 
papers and from our study. Our results are labeled with 
asterisks. So, in this figure, we have seventeen results 
for the flux density values S v of Cygnus loop derived in 
several papers (values are given in Table 3 of Uyaniker 
et al. [ll| and five results from our paper. 

It can be seen from Fig. [IT] that our added flux values 
fit very well among other fluxes, which shows correctness 
of our method, and also that we determined loop area 
and its brightness well. 



TABLE V: Flux densities (Jy) that we calculated at the five 
frequencies 



Frequency (MHz) 


Flux density (Jy) 


2720 


83.37 ± 2.60 


1420 


80.05 ± 8.13 


820 


139.35 ± 12.10 


408 


235.44 ± 15.36 


34.5 


1309.04 ± 65.65 



V. CONCLUSIONS 

The main result is method for determination bright- 
ness temperature given in article Borka [25| which is 
developed for large radio loops, and we show that this 
method is very good for much smaller loop, e.g. Cygnus 
loop. We check the method by applying it to the Cygnus 
radio loop at the wide range of frequencies. It is in good 
agrement with method of T — T graphics and results 
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FIG. 11: Spectrum of Cygnus Loop: flux versus frequency 
for more different frequencies, obtained from all available flux 
density values (values are given in Table 3 in Uyaniker et al. 
(2004)) and from values calculated in this paper. We identify 
our values by asterisks. 



obtained with another methods. This method is quite 
simple because we use brightness temperature isoline to 
define border of a Cygnus loop. Other authors are using 
different squared or rectangular areas to determine area 
of the loop and calculate spectral indices, brightness tem- 
perature and the flux density of the loop (Uyaniker et al. 
[TTI | . page 917 and Leahy & Roger [nj], page 786). Also, 
we calculated flux for Cygnus loop and compare our re- 
sults with results of other authors in different ranges of 
frequencies. We show that our results are in good agre- 
ment with these results. 

We estimated the temperatures and brightness of the 
Cygnus loop SNR on the basis of observations of the con- 
tinuum radio emission at the frequencies: 2720, 1420, 
820, 408 and 34.5 MHz (this paper and Borka Jovanovic 
& Urosevic [22[ ). The sensitivity of the brightness tem- 
peratures are: 5 mK for 2720 MHz, 50 mK for 1420 MHz, 
0.2 K for 820 MHz, 1.0 K for 408 MHz and about 700 
K T b for 34.5 MHz. At the frequency of 2720 MHz the 



measurements are the most precise (with the least rel- 
ative errors) so positions of the brightness temperature 
contours of the loop are the most realistic for this fre- 
quency. 

Borders between the five frequencies are somewhat dif- 
ferent for this SNR probably due to small, random and 
systematic errors in the calibrated data. Also, we sup- 
pose there are uncertainties of about (2 x AT) 10 mK for 
2720 MHz, 100 mK for 1420 MHz, 0.4 K for 820 MHz, 
2.0 K for 408 MHz and about 1400 K T b for 34.5 MHz, 
in the border due to measurement errors, and there is a 
tiny difference in the absorption of radio emission in the 
interstellar medium at different wavelengths (Pacholczyk 

We determined average brightness temperature from 
a Cygnus radio loop region, after subtraction of a back- 
ground level. Our obtained values (when all reduced to 
1000 MHz for comparison) are in good agreement with 
the earlier results. We present the radio continuum spec- 
trum of the Cygnus loop using average brightness tem- 
peratures at five different frequencies. As it can be seen 
from Figs. [7]and[51 given linear fit provides reliable spec- 
tral index. We present the T — T plots which enable the 
calculation of spectral index, too. 

Also, from our results can be concluded that Cygnus 
loop evolves in the low density environment and the ini- 
tial energy of SN explosion was relatively low. This can 
be read after superposing the position of this loop to 
theoretical E — D diagram form Berezhko & Volk |38| . 
which was derived from non-linear kinetic theory of dif- 
fuse shock acceleration mechanism. 

We showed that method for defining a loop border and 
for determining the values of temperature and brightness, 
which we developed for main Galactic loops I- VI, could 
be applicable to all SNRs. 
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